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ABSTRACT 

In recent years increasing evidence has emerged for a thermal component in the 7- and X-ray 
spectrum of the prompt emission phase in gamma-ray bursts. The temperature and flux of the 
thermal component show a characteristic break in the temporal behavior after a few seconds. We 
show here, that measurements of the temperature and flux of the thermal component at early times 
(before the break) allow the determination of the values of two of the least restricted fireball model 
parameters: the size at the base of the flow and the outflow bulk Lorentz factor. Relying on the 
thermal emission component only, this measurement is insensitive to the inherent uncertainties of 
previous estimates of the bulk motion Lorentz factor. We give specific examples of the use of this 
method: for GRB970828 at redshift z — 0.9578, we show that the physical size at the base of the flow 

is r = (2.9 ± 1.8) x 10 s Y~ 3/2 cm and the Lorentz factor of the flow is T = (305 ± 28) r o 1/4 , and for 

GRB990510 at z = 1.619, r = (1.7 ± 1.7) x 10 s Y~ 3/2 cm and T = (384 ± 71) Y Q 1/4 , where Y = 1Y 
is the ratio between the total fireball energy and the energy emitted in 7-rays. 

Subject headings: gamma rays: bursts — gamma rays: theory — plasmas — radiation mechanisms: 
non-thermal — radiation mechanisms: thermal 



1. INTRODUCTION 

In recent years increasing evidence has appeared that, 
during the first stages of the prompt emission of long 
duration gamma-ray bursts (GRBs), a thermal compo- 
nent accompanies th e underlying non- t herma l emission 
iRvde I (12004 120051); [gampana et ~aT\ ( 1 2Q06J)- see also 
iGhirlanda et al. I (|2003t ): iKaneko et al. I (120031 )]. 5 An 
analysis of BATSE bursts that are do minated by quasi- 
thermal emission (jRvde 1 12004 |2005[ ) showed that the 
observed temperature exhibits a similar behavior in all 
of them: an initially (approximately) constant temper- 
ature at a canonical value Tq 6 - ~ 100 keV which af- 
ter ~ 1 — 3 seconds decreases as a power law in time 
T ob - cx t~ a , with power law index a ~ 0.6 — 1.1. The 
redshifts of most of these bursts are u nknown. An ad- 
ditional analysis (|Rvde fe Pe'er I [20071 ) shows that after 
a short rise, the flux of the black body component of 
these bursts also decreases with time as Fg b B cx t~@ , with 
(3 ~ 2.0 — 2.5. We showed there that this temporal behav- 
ior can be ex plained as due to the h i gh latitude emission 
phenomen on (|Fenimore et al. I fl996: Gra not et al. Ill999l : 
[Qm1l2002l) . 

According to the standard fireball scenario, the non- 
thermal photons originate from the dissipation of the fire- 
ball kinetic en ergy. The dissipation mechanism (e.g., in- 
terna l shocks (|Paczvnski fc Xulll994 iRees fc Meszaros 1 
|1994T >. magnetic reconnection ([Giannios fc Spruit 1 120051 : 
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Gianniosl I2006D or external sh ocks (Meszar os fc Reesl 
19931: ID ermer fc Mitmanl ll999)) is yet uncertain, and 
can in principle occur at various locations. As opposed 
to this ambiguity in understanding the origin of the non- 
thermal component, the thermal component must origi- 
nate at the photosphere. According to the high latitude 
emission interpretation of the data, the highest temper- 
ature and the maximal thermal flux initially observed 
are emitted from the photosphere on the radial axis to- 
wards the observer. Thus, in principle the radius of the 
emission site of these photons can be determined. 

In this Letter we show that combined early time mea- 
surements of the observed temperature and thermal flux 
for bursts with known redshift allow us to directly deter- 
mine the values of the bulk motion Lorentz factor, the 
physical size at the base of the flow and the photosphcric 
radius. This is due to the fact that the observed tem- 
perature and flux of the thermal component depend on 
three internal parameters only: the isotropic equivalent 
luminosity of the thermal component Lbb, the Lorentz 
factor of the bulk motion of the flow at the photosphcric 
radius r\ and the physical size at the base of the flow ro, 
and that Lbb can be directly measured for bursts with 
known redshift and measured thermal flux. In fj2] we 
give a short description of the model, and implications 
are given in In 21 we summarize and compare our 
results with those of previous methods of estimations of 
the bulk motion Lorentz factor. 

2. MODEL: EXTENDED PHOTOSPHERIC EMISSION 

In the cl a ssical fireball model of g amma-ray bursts 
(iGoodman I 119861 : IPaczvnski I Il986t IPaczvnski fc Xu I 
|1994|) . a thermal plasma of electrons, positrons, and pho- 
tons expands rapidly from an initial radius ro . Conserva- 
tion of energy and entropy imply that the bulk Lorentz 
factor of the flow increases as Y(r) cx r, until the plasma 
reaches the saturation radius r s — tjtq, above which the 
plasma Lorentz factor coasts with T = 77 = L/Mc 2 . 
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Here, L is the isotropic equivalent burst luminosity, M 
is the mass ejection rate and c is the speed of light (from 
here on we restrict the discussion to long bursts, charac- 
terized by extended emission of relativistic wind). 

The photospheric radius r p h is the radius above which 
the flow becomes optically thin to scattering by the 
baryon related electrons. Depending on the values of 
the free model parameters (77 , L and rp) this radiu s can 
be smaller or larger than r s (Meszaros et al. 2002): for 
77 > (<)t7* = (Lcr T /47rm p c 3 r ) 1/4 , r ph < (>)r s . Here, a T 
is the Thomson cross section and m p is the proton mass. 
The luminosity L is measured for bursts with known red- 
shift, L = And 2 L Y F ob - , where d,L is the luminosity dis- 
tance, F ob - is the total (thermal + non thermal) observed 
7-ray flux, and Y = e/e 7 > 1 is the ratio between the 
total fireball energy and the energy emitted in 7-rays. 
As we show below, the measurement of ro is similar in 
both scenarios, r p h > (<)r s . Thus, it is possible to deter- 
mine whether r p h is below or above r s from measurable 
quantities, and to determine 77 in the second case. 

The thermal component originates from the photo- 
sphere of an expanding plasma jet. The observed ther- 
mal flux (integrated over all frequencies) is given by in- 
tegrating the intensity over the emitting surface, F^g — 
(27r/d|) / d/j,fxrl h V 4 (aT' i /'!r). Here, V is the comov- 
ing temperature at the photospheric radius, a is Ste- 
fan's constant and T> = T>(0) is the Doppler factor, 
V = [r(l — /3/j)]" 1 . The angle 9 is the angle between 
the direction of the outflow velocity vector (/3) and the 
line of sight, (i = cos(6), and T = (l-/? 2 ) -1 / 2 is the out- 
flow Lorentz factor. For a plasma Lorentz factor much 
larger than the inverse of the jet opening angle T 3> 9J 1 
and for early enough times, at a given observed time the 
integration boundaries are determined uniquely by the 
emission duration, regardless of the value of 9j. 

Due to the Doppler effect and to the cosmological red- 
shift, photons that are emitted at frequency v 1 in the co- 
moving frame of a relativistically expanding plasma with 
Lorentz factor T at redshift z, are observed at frequency 
v ob - = V(6)v' /(l + z). For T > 1 and emission on the 
line of sight, v ob - ~ 2lV/(l + z). In the exten ded emis- 
sion interpretation of the data (|Rvde fe Pe'er II2007T ). an 
observer sees simultaneously photons that originate from 
a range of angles to the line of sight, 9 min < 9 < # max . 
In this case, a thermal spectrum (with temperature T 1 ) 
in the comoving frame is observed as a modified black- 
body spectrum. Nonetheless, the observe d spectrum i s 
very similar to a pure black-body spectrum (jPe'er Il20"07f ). 

During the first few seconds, when the observed tem- 
perature is nearly constant, the observed radiation is 
dominated by photons emitted close to the line of sight, 
i.e., # m in = 0. According to this interpretation, at 
t ob ' > tbreak-, there is no more emission from 9 = 
because the inner engine activity decreases, and the 
emission is dominated by high latitude effects. Here, 
tbreak is the break time in the temperature's temporal 
behavior. At t ob - < i&reafc, the observed spectrum is 
very close to a black bo dy spectrum with temperature 
T ob - ~ 1.48IT7(1 + z) (|Pe'er 112007ft 6 . During this pc- 

6 The relation between T and T' is often written in the lit- 
erature as T ob - ~ TT'/(1 + z). In fact, T>{0 = 1) ~ 2T, thus one 
should add an extra factor of 2. The factor 1.48 used here results 



riod, the upper integration boundary in the equation for 
the thermal flux is /i max = cos(# m i n ) = 1, and the ratio 

(- F Bs/ crT ° b ' 4 ) 1/2 which we denote as TZ, is equal to 
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The prefactor (1.06) originates from the dependence of 
the photosph eric radius on the angle to the line of sight 
(jPe'er 1120071) . 

We can now make the discrimination between the two 
possible cases: r p h < (>)r s . If r p h < r s , then Y{r) oc r. 
In this case, r p h/T = ro, and equation [T] becomes 

ro ^ h<r ' )= (im) {i + zf n - (2) 

In this case, it is not possible to determine the photo- 
spheric radius, or the value of T(r p h). 

In the case r p h > r s , the photospheric radius is 
given by r vh = (La T /8irr] 3 m v c 3 ) ()Meszaros et al. Il2002l : 
iDaigne fc Mochkovitchll2002t iBroderick ll2005t )~At this 
radius, the comoving temperature is given by 

1/4 / \ -2/3 
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where a is the radiation constant 

7>ob. 



Setting r = 77 and L = 
ATrd 2 L Y F ob - in the equation of the photospheric radius, 
one obtains from equation Q] the coasting value of the 
Lorentz factor, 



77 = 



(1.06)(1 



2 YF° b a T 
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ZrripC 6 K 
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Equations Q] IS1 and [H now give the physical size at the 
base of the flow, 
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Mr P H > r.) = n J 6n nR , in ^ ( ) n . 

(5) 

We thus find that a measurement of TZ and the ratio 
of the black body flux to the total flux at the very early 
observed times from bursts with known redshift give a 
direct measurement of ro, and that the result is similar 
(up to a numerical factor of the order unity, provided that 

Y is not much larger than 1; see discussion on the value of 

Y in Q, for the two considered cases, r p h < (>)r s . The 
measured values of ro and L can be used to determine the 
value of 77*, which is independent of the specific scenario. 
One can then use the measured values of TZ and F ob - to 
determine the value of 77 using equation^ If the obtained 
value is larger than the value of //» , then r p h < r s , in 
which case equation |4] should not be used and the value 
of 77 remains undetermined. 

3. IMPLICATIONS 

Relations |4] and [5] allow a direct measurement of the 
size at the base of the flow, ro and of the bulk motion 
Lorentz factor of the flow from GRBs with known red- 
shift and energy content. In addition, it is possible to 
determine the photospheric radius r p h and the satura- 
tion radius r s , if measurements of the thermal flux and 
temperature are available at early enough times. 

from the angular integration. See lPe'eTI 1 120071) for details. 



(1.48) 6 (1. 06) 4 (1 + z) 2 \YF ob 
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Parameters determination using thermal component 
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We illustrate the use of this method on two bursts with 
known redshifts, namely GRB970828 and GRB990510 
observed by the BATSE detectors aboard the Compton 
Gamma Ray Observatory. BATSE detected bursts in the 
20 keV - 2 MeV energy range. The time-resolved spectra 
for the selected bursts were fitted using a Planck function 
and a single power-law to model the photospheric and the 
non-thermal emission com ponents, following the metho d 
presented in lRvde I ([20041) (see also lRvde fe Pe'er Il2007t ). 

The analysis of the thermal component of GRB970828 
is presented in figures [TJ This burst, at redshift z = 
0.9578, had a good temporal coverage for the first 100 
seconds. The left hand panel in Figure Q] shows the tem- 
poral behavior of the temperature of the thermal com- 
ponent. During the first ~ 8 seconds the observed tem- 
perature rises slightly to a value of 78.5 keV, after which 
it shows a rapid decrease that can be fitted as a power 
law in time with a power law index a = —0.51. In the 
right hand panel of Figure [T] we show the temporal be- 
havior of the function 1Z. This function shows a slight 
increase during the first ~ 7 seconds, after which it rises 
as a power law in time with a power law index f3 = 0.67. 
The smooth increase of 1Z before the break implies that 
there is no significant energy dissipation below the pho- 
tosphere, which, if it occurred would result in a strong 
fluctuation and affect the smoothness of 1Z. The break 
times in the temporal behavior of the observed temper- 
ature and 1Z are the same within the errors. 

According to the extended high latitude emission inter- 
pretation of this result, we deduce that the first episode 
of significant inner engine activity took place during the 
first 7 — 8 seconds, and at later times we are observ- 
ing photons emitted off axis (we neglect here late time 
episodes of engine activity that occur after ~ 25 s and 
~ 60 s in this burst). The values of ro and 77 are cal- 
culated using the observed values of the temperature 
T ob - = 78.5 ± 4.0 keV , 1Z = (1.88 ± 0.28) x 10~ 19 , and 
the ratio of thermal to total flux F° b B /F° b - = 0.64 ±0.20 
at the break time. The error bars on the measured 
quantities are averaged over the first seconds, before 
the temporal break. Considering a flat universe with 
SIa = 0.73, Hq = 71 km/s/Mpc, the luminosity distance 
for this burst is — 1.94 x 10 28 cm. Using equa- 
tions H and we find that T = (305 ± 28) F 1/4 and 
r Q = (2.9 ±1.8) x 10 8 Y~ 3/2 cm, where Y = 1Y . The cal- 
culated value of rj* = 463 Y^ 8 proves that indeed r p h = 



2.7 x 10 11 F 1/4 cm is larger than r s = 9.0 x 10 lu Y^ '* cm, 

which implies that r/ = T = 305 Yq' 4, is the coasting value 
of the outflow Lorentz factor. The statistical error on the 
estimated value of 77 is < 10%. The systematical error 
results from the uncertainty in the value of Y , and is not 
expected to be more than tens of percents (see $4] be- 
low), giving the best constraint on the estimated value 
of rj measured so far. 

A similar analysis was carried out for GRB990510 at 
z = 1.619. The results obtained are similar, ro = 

(1.7 ±1.7) x 10 8 V 3/2 cm and X] = (384 ± 71) Y" 1/4 . The 
larger statistical errors compared to GRB970828 mainly 
reflect the fewer available data points for this burst. 

The value of 77* = 830 Y^ 8 proves that indeed r p h = 

7.7 x 10 11 F 1/4 cm is larger than r s = 6.3 x 10 10 F~ 5/4 . 
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Fig. 1. — Temporal behavior of the thermal component in 
GRB970828 at 2 = 0.9578. Left panel: the observed temperature. 
During the first ~ 8 s the temperature rises slowly up to 78.5 keV, 
after which it decreases as a power law in time with power law 
index —0.51. The fit was made on the data up to 20 s, which in- 
cludes the first pulse structure of the light curve. Right panel: the 
temporal behavior of the ratio Tl = {F BB / ^T obi ) 1 / 2 . This ratio 
increases slowly during the first ~ 7 s, after which it increases as 
a power law in time with power law index 0.67. The break time 
after ~ 7 s is close to the break time in the temporal behavior of 
the temperature. 

The temporal behavior of the observed temperature and 
TZ were found to be similar in a large sample of BATSE 
bursts, providing further evidence for our model. How- 
ever, the redshift of most of these bursts is unknown, 
thus definite values of 77 a nd r n could not be ob tained. 
The full sample appears in lRvde fc Pe'er I (|2007fh 

4. DISCUSSION 

In this Letter we showed that by measuring the ob- 
served temperature and thermal flux of the thermal com- 
ponent that accompanies the prompt emission of GRBs, 
it is possible to determine the values of two of the least 
restricted parameters of the fireball model: the size at 
the base of the flow and the outflow bulk Lorentz factor 
(in the case that the photospheric radius is larger than 
the saturation radius). In this case, it is also possible to 
determine the saturation radius r s and the photospheric 
radius r p h- We showed that the calculation of the ini- 
tial size of the flow (equations [2] [5]) is similar (up to a 
constant of the order unity) for the cases r p h > (<)r s . 
This allows a comparison between the measured value of 
r\ (equation 2]) and the derived value of 77*, and a discrim- 
ination between the two cases. We have given examples 
of the use of this method for the determination of r\ and 
ro for two specific GRBs. 

The largest uncertainty in the estimate of 77 is due 
to the uncertainty in the value Y. Value of 1 < 

Y < 3 — 5 were suggested based on afterglow obser- 
vations dFreedman fc Waxman I [2001b iFrail et al. I [2001b 
iFriedman fc Bloom 1 120051 : iGranot et al. 112006ft Theo- 
rctical arguments based on fitti ng the flux of ultra -high 
energy cosmic rays (UHECR; IWick et al~\ l2004ft . un- 
der the assumpt i on that UHE CRs originate from GRBs 
(jWaxman 1995: Vietri 199 5]), a s well as efficiency con- 
siderations (|Fan fc Piran I [2006) suggest larger value, 

Y > 10. For bursts with a dominant thermal com- 
ponent, as considered here, we expect Y to be close 
to unity. If afterglow measurements of bursts with de- 
tected thermal component and known redshift become 
available, as expected after the launch of the GLAST 
satellite, this uncertainty could be removed. A second 
source of systematic uncertainty can result from domi- 
nating Compton scattering (which conserves the number 
of photons) resulting in a Wien spectrum, rather than a 
thermal spectrum. Observationally, a Wien spectrum is 
hard to discriminate from a thermal spectrum. In this 
case, the systematic error in estimating the temperature 
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is (3k B T)/(2.7k B T) ~ 10%, which transforms into - 5% 
uncertainty in the estimated value of r/. Currently, this 
uncertainty is smaller than the statistical uncertainty, 
> 10%. 

Other methods of estimating the bulk motion Lorentz 
factor in GRBs relied on a large number of uncertain 
model assumptions and uncertainties in the values of the 
free model parameters. A widely used lower limit for 
77 is obtained by calculating the minimum Lorentz fac- 
tor required in or der for the observed energetic photons 
not t o annihilate (iKrolik fe Pier 1 119911: iFenimore et al.\ 



1993; Wo ods fc Loebl 119951 : iBaring fc Harding I 119971 : 
Lithwick fc Sari 1120011 ). In addition to providing only a 
lower limit, in order to get a good estimate of 77 m i n a wide 
spectral coverage of the GRB emission, from the optical 
band to the 7-rays is required. In these calculations, the 
prompt emission spectru m is sometimes approxim ated as 
a broken power law (e.g.. iLithwick fc Sari 1120011). w hich 
may be too simplified (e.g. JPe'er fc Waxman MM)- 

An alternative method to estimate the Lorentz factor 
is by modeling the early afterglow emission, on the as- 
sumption that the optical flash observed in a few cases 
results from synchrotro n emission by electr ons heated 
by the reverse shock (|Sari fc Piran ] H999I) . A seri- 
ous drawback of this method is that the estimate re- 
lies on the poorly known shock microphysics parame- 
ters (such as e P ., en etc.). A mor e advanced method, 
introduced by IZhang et al. I (|2003l ). relies on compar- 
ing the emission from the forward and reverse shock 
during the early afterglow. An underlying assump- 
tion in this estimate is that the values of the micro- 
physics parameters at the forward and the reverse shocks 
are similar. Other methods rely on measurements of 
the physical parameters during the late afterglow emis- 
sion, assuming that the flow expands in a self-similar 
motion during this phase. The initial value of the 
Lorentz factor is deduc ed by measuring the rise time 
of the early afterglow (ISari 1 119991: IWang et al. 1 120001: 



Soderbcrg & Ramirez-ruiz I 120021 iKobavashi fc Zhang I 
2003). An inherent drawback of this method is the as- 
sumption that the microphysical parameters are constant 
in time during the late afterglow. 
The method presented here of estimating r\ is indepen- 



dent of any of the uncertainties inherent in the former 
methods. Moreover, it gives a direct measurement of rj, 

rather t han a lower lim it. The results presented in §3\ 

(see also lRvde fc Pe'er~ll2007D indicate values of 77 close 
to the earlier estimates. In addition, the values found 
for the size at the base of the flow ro could further con- 
strain GRB progenitor models. The statistical errors on 
the values of these numbers are much smaller than any 
previous estimates. 

These facts have several important consequences. 
First, they strengthen the interpretation of the prompt 
emission as being composed of a thermal component, in 
addition to the non-thermal component. Therefore, any 
interpretation of the prompt emission data must take this 
thermal component into account. Second, it shows that 
the extended high latitude emission interpretation of the 
late time temporal behavior of the thermal component is 
consistent with the fireball model predictions. Thus, this 
interpretation may also be used to understand the strong 
X-ray flares observed by the SWIFT satellite. Third, the 
consistency found between the different methods for es- 
timating the value of rj can be used to strengthen the 
validity of the underlying assumptions in previous esti- 
mates of i], such that the values of the microphysical pa- 
rameters (e e , 6b, etc.) are indeed constant in time during 
the afterglow emission phase. And last, the direct mea- 
surement of the physical size at the base of the flow is 
another, independent indication that a massive star is 
indeed the progenitor of long duration GRBs. Since rj 
is related to the mass ejection rate, our measurements 
could be useful to constrain models of GRB progenitors. 
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